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Abstract 

O t 

We investigate magnetic properties of color superconducting quark matter 
within a Ginzburg-Landau approach. The simultaneous coupling of the quark 
fields to gluonic and electromagnetic gauge fields leads to rotated electromag- 
netism with a massive (Higgsed) and a massless photon-gluon field. We derive 
the Ginzburg-Landau equations for superconducting quark matter taking into 
account the rotated electromagnetism in the general case when the rotation 
angle a is an arbitrary function of the coordinates. We solve these equations 
for the isolated vortex in superconducting quark matter. We obtain a solution 
for the magnetic and gluomagnetic fields and expressions for the calculation 
of the penetration depth and the quantized magnetic flux in quark matter. 
For this case we have demonstrated that the occurence of electric and color 
Meissner currents is a consequence of the color superconducting state of quark 
matter. 
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I. INTRODUCTION 



Recently, the possible formation of diquark condensates in QCD at finite density has 
been reinvestigated in a series of papers following Refs. |l]J3]. It has been shown that in 
chiral quark models with a nonperturbative 4-point interaction motivated from instantons 
or nonperturbative gluon propagators |4],|| the anomalous quark pair amplitudes in the 
color antitriplet channel can be very large: of the order of 100 MeV. Therefore, in two-flavor 
QCD, one expects this diquark condensate to dominate the physics at densities beyond the 
deconffnement /chiral restoration transition and below the critical temperature 50 MeV) 
for the occurence of this two-flavor color superconductivity (2SC) phase. In a three-flavor 
theory it has been found |||7j that there can exist a color-flavor locked (CFL) phase for not 
too large strange quark masses M where color superconductivity is complete in the sense 
that diquark condensation produces a gap for quarks of all three colors and flavors, which 
is of the same order of magnitude as that in the two-flavor case. 

The high-density phases of QCD at low temperatures are most relevant for the explana- 
tion of phenomena in rotating compact stars - pulsars. Conversely, the physical properties 
of these objects (as far as they are measured) could constrain our hypotheses about the 
state of matter at the extremes of densities. In contrast to the situation for the cooling 
behaviour of compact stars where the CFL phase is dramatically different from the 2SC 
phase 0, we don't expect qualitative changes of the magnetic field structure between these 
two phases. Consequently, we will restrict ourselves here to the discussion of the simpler 
two-flavor theory first. The 2SC phase window is expected to occur at densities just above 
the deconfinement transition inside a compact star where the strange quark flavor is either 
still confined [10| or the self consistently determined strange quark masses are large enough 
to entail color-flavor unlocking . 

According to Bailin and Love Jl2| the magnetic field of pulsars should be expelled from 
the superconducting interior of the star due to the Meissner effect and decay subsequently 
within « 10 4 years. For their estimate, they used a perturbative gluon propagator which 
yielded a very small pairing gap and they made the assumption of a homogeneous magnetic 
field. Since both assumptions seem not to be valid in general, we have performed in Ref. 
|T3| a reinvestigation of the question with the result that 2SC quark matter is a type-II 
superconductor which can form vortices in response to the external magnetic field thus 
behaving similar to protons in the superfluid neutron phase of a neutron star |14,1^|. The 



magnetic field in the quark core of the neutron star can then exist for periods much longer 
than the spin-down age of a pulsar. 

The authors of Ref. [16| have considered this question taking into account the rotated 
electromagnetism. They came to the conclusion that the magnetic field will live in the quark 
core sufficiently long although it does not form the structure of quantum vortices because it 
obeys the force-free Maxwell equation. 

In the present work we derive the Ginzburg-Landau equations with account of rotated 
electromagnetism and solve these equations for the example of an isolated quantum vortex 
in quark matter which demonstrates that there is always an electric and color Meissner 
current when the quark matter is superconducting. This conclusion holds also for the case 



of a homogeneous external magnetic field as we have shown in a separate paper JTTj 
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II. GINZBURG-LANDAU EQUATIONS 



In the paper [O] we have obtained the Ginzburg-Landau equations for relativistic super- 
conducting quarks, supposing that the superconducting quark matter phase with ud diquark 
pairing has P p = 0, where p is the color antitriplet index. We study the consequences of 
superconducting quark cores in neutron stars for the magnetic field of pulsars. We find 
that within recent nonperturbative approach to the effective quark interaction the diquark 
condensate forms a type II superconductor whereas previously quark matter was consid- 
ered as a type I superconductor [|12[]. In both cases the magnetic field which is generated 
in the superconducting hadronic shell of superffuid neutrons and superconducting protons 
can penetrate into the quark matter core since it is concentrated in proton vortex clusters 



where the field strength exceeds the critical value. Recently, in the paper [16[], discussing 
the Meissner effect for color superconducting quark matter, the authors introduce the Q- 
charge generator 

Q = Q + vPs, (1) 

where Q is the conventional electromagnetic charge generator and Pg is associated with the 
one of the gluons in the representation of the quarks 

Q = diag (—, ——, — -\ in flavor u, d, s space , (2) 

\ O O O / 

Pg = -^diag (1, 1, —2) in color r, g, b space . (3) 
V3 

The condition they requested was that the Q- charge of all Cooper pairs which form the 
condensate vanishes 

Qd p = . (4) 

Applying this condition (instead of Pg = 0) for the derivation of the Ginzburg-Landau 
equations, we will see that in the superconducting quark matter two different types of 
magnetic fields can exist simultaneously. The penetration depths X q for these two magnetic 
fields are completely different. For one of them the penetration depth is finite and much 
shorter than has been found in |13|], but for the other one it becomes infinitely large. 



The free energy density in the superconducting quark matter with ud diquark pairing 
can be written in the following form |J 



/ = /„, + ad p d* p + b(d p d;f + 7 |Prf p | 2 + + , (5) 



where a = dn/dE t, (3 = dn/dE 7((3)(nk B T c )- 2 /8, 7 = p| /?/(6/i 2 ) and t = (T - T c )/T c 
with T c being the critical temperature, pp the quark Fermi momentum. In zeroth order with 
respect to the coupling constant, dn/dE = hpf/tt 2 . \Pd p \ 2 is the kinetic term of the free 
energy density (|5|) which is given by 

\Pd p \ 2 = |(V + xeAQ + igG 8 P 8 )d p \ 2 . (6) 
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If we take into account the condition (Q) and use (|l|) we will have 



Qd = - V P 8 d . (7) 

Inserting (^) into (|6|) we will get 

| Ay 2 = |(V + ierjAP 8 + tgG 8 P s )d p \ 2 . (8) 



Following the paper |lq| , let us introduce instead of the original gauge fields A and G 8 the 
linear combinations A x and A y 

A x = — =====- = - sin aA + cos aG 8 , (9) 
Vr/^e^ + g z 

X gA + rjeGs - . - . , 

A„ = —===== = cos a A + smaGs . (10) 

V^ 2 e 2 + g 2 

This corresponds to a rotation with the angle a in the orthogonal basis of the fields A x and 
A y , where 

while the normalization is preserved 

Al + Al = A 2 + Gl. (12) 

At neutron star densities the gluons are strongly coupled (g 2 /4ir w 1) and the photons 
are of course weakly coupled (e 2 /47r w 1/137), so that a w ^e/^ is small. For the diquark 
condensate, where blue-green and green-blue ud quarks are paired, r\ — 1/ \f?> and therefore 
a ~ 1/20. Taking into account © and (|T0|) , the kinetic term (§) will become 

\Pd p \ 2 = (V - iqT]A x )d* p (V + iqr]A x )d p , (13) 

where 



q=y/ri*e* + g*P 8 , (14) 

and P 8 = l/v^- We see that the kinetic term has the conventional expression but instead 
of the charge of the ud diquark pair equal to e/3 we have the new charge q which is much 
larger than e/3 (about 20 times). The Ginzburg-Landau equations are obtained in the usual 
way. If we demand that the variation of the free energy with respect to the parameters d*, 
A x and A y have to be equal to zero, we will get 

= ad p + p(d p d* p )d p + 7(iV + qA x fd p , (15) 
the equation of motion for the diquark condensate and for the gauge fields 
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= sin a rot rot A — cos a rot rotGg — 4iriq'y[dVd* — cTV<i] — 8nq 2r y\d\ 2 A x , (16) 
= cos a rot rot A + sin a rot rotGg • (17) 

These equations can be written in the following form 

x2 + .7 . • 2 7 . sina {d p Vd* - d*Vd p ) . 

A,, rot rot/L + sin aA = i— h sinacosaGs (18) 

y 2g \d\ 2 

\2 + , 2/1 .cosa(d p Vd* - d*Vd p ) . 

A„ rot rotG 8 + cos aG 8 = — z h sin a cos ct/L , (19) 

y 2q \d\ z 



where 

V = y[toH\d\q 



Aq 2 tp 3 F 



(20) 



This is the system of Ginzburg-Landau equations for superconducting quark matter which 
takes into account the "rotated electromagnetism" . When we define in the equations fll8|) 
and ( |r9| ) the order parameter in the following form 

d p = \dp\e** , (21) 

where <fi is the phase of the order parameter, then this system of equations can be written 
in the following form 

— * — * $ Sill Oi —* 

A 2 rot rotA + sin 2 aA = i— V0 + sin a cos aG 8 (22) 

2n 

A 2 rot rotG 8 + cos 2 aG$ = —i— V0 + sin a cos aA (23) 

2n 

Here $ g has the well known form 

2nhc 

$ 9 = — - . (24) 

If we take the contour integral on both sides of equation (|22"D where the contour L is chosen 
as the equatorial ring which limits the quark matter core region and take into account that 

rot rot! dl= I rotB dl = , (25) 

L J L 

G 8 dl = , (26) 

L 

so that the electrical currents and the vector potential of the gluon field on the surface of 
the quark core must vanish, then we get 

A dl= [ rot A dS= [ B dS = -^-N . (27) 
l Js Js sin a 
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This means that the flux of the magnetic field through an arbitrary surface over the contour 
L is equal to the sum of fluxes which are generated by N vortices each of which has a flux 
$ 9 /sina . 

When we calculate this flux we obtain 

, 6irhc , , 

= — = 6 $ , (28) 

where $ = 2 x 10~ 7 G cm 2 is the magnetic flux quantum. In order to obtain this result 
we have used equation fl22|). It is easy to see that the same result can be obtained using 
equation (|23|), as to be expected. 



III. SOLUTION OF GINZBURG-LANDAU EQUATIONS FOR ISOLATED 

VORTEX IN QUARK MATTER 

We now consider the solution of the Ginzburg-Landau equations |23| ) for an isolated 
vortex situated in homogeneous superconducting quark matter. To this end it is sufficient 
to find the solution of these equations without the vortex terms. The contribution of an 
isolated vortex is taken into account by the condition ( ^7|) where N = 1. Consequently we 
have to solve the following system of equations 

A 2 rot rot A + sin 2 a A = sin a cos a G$ (29) 
A 2 rot rotG 8 + cos 2 a G$ = sin a cos a A . (30) 

From these equations it is easy to see that 

rot rotGs = — cot a rot rotA . (31) 

When we act on both sides of equation (|29"D with the operator rot rot and take into account 
that a = const, we obtain 

A 2 , rot rot rot lotA + sin 2 a rot rotA = sin a cos a rot rotGs • (32) 

Inserting (|3lD into (|32|) , we obtain the equation for the determination of the electromagnetic 
vector potential 

A 2 rot rot rot rot A + rot rot A = . (33) 

Let us denote 

rot rot A = M . (34) 

Then, equation (31) takes the form 

A 2 , rot rotM + M = . (35) 
Consequently, instead of equation ( |33"D we can solve the system of equations ([!4|) and (|35|). 
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Let us consider the solution of equation ( |35|) . Since our problem has cylindric symmetry, 
the unknown function shall depend only on the coordinate r denoting the distance to the 
vortex. It is easy to see that the vectors M and A have only azumutal components M^(r) 
and A v (r). Then equation (^) takes the form 

^ + ^-<H )M ^°' (36) 

The general solution of this equation which tends to zero at infinity is given by 

M ip (r) = cK 1 (r/\ q ) , (37) 



where K\ is the modified Bessel function of first kind. Next we solve the equation fl34|) . It 
is easy to see that a special solution of this equation is A = —X q M, and the general solution 
which vanishes at r — ► oo will be C2/r . Finally for the vector potential we obtain 

A v {r) = c x K 1 {r/\ q ) + c 2 /r . (38) 

The nonvanishing component of the magnetic field has z- direction and is determined by the 
formula 

B t (r) = ~(rMr)) ■ ( 39 ) 
r ar 

Inserting fl3~8"D into fl39|) we finally obtain 

B g (r) = a K (r/X q ) , (40) 

where Kq is the modified Bessel function of zero kind. The constant C\ is determined from 
the equation (|27j) as 

$' 

- 7db ■ ( 41 ) 



Then we have 



r 



b -<- t) = *m k ° r • (42) 
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If we introduce the gluomagnetic field B 8 = rotG$ , then we can determine it using equation 
(pl|). The condition that the magnetic fields B and B& have to vanish at r — > cxo, immediately 
determine the solution for B$ 

B 8 = -B cot a, (43) 

where 5 is given by equation fl42f) . 

As can be seen from the solution (E2|) the field of the quantum vortex in quark matter 
formally has the same form as the field of a vortex in an ordinary superconductor. Whereas 
in the case of the ordinary superconductor in the definition of the flux $ and the penetration 
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depth A enters the same charge of the Cooper pair (i.e. 2e, where e is the charge of the 
electron), in quark matter in the definition of $ g enters the charge of the quark Cooper 
pair e/3, and therefore $' q = 6$ , an d in the definition of the penetration depth X q enters 
the charge q, which is about 20 times larger than e/3. Therefore, the penetration depth X q 
which is proportional to q~ x is 20 times smaller than that of hadronic matter. It is obvious 
from equation (^) that the magnetic quantum vortex is accompanied by a gluomagnetic 
quantum vortex the center of which coincides with that of the former, but its force lines 
have the opposite direction to the magnetic ones. The flux of the gluomagnetic vortex is 
greater than that of the magnetic vortex by a factor cot a. 

Let us consider the equations for the rotated magnetic fields B x = rotA x and B y = rotA y . 
which we obtain from the equations (|i"6|), ( |17|) 

X 2 q rot rot4 + B X = %e5(r) , (44) 
rotBy = . (45) 

The energy of isolated quark vortex is given by 

E = ±- f[B 2 y + B 2 x + (X q rotB x ) 2 ] dV (46) 

o7T J 

Note that from equations fllOD and (f43|) follows that the magnetic field B y of the isolated 



quark vortex is equal to zero. The magnetic field B x will be found from the solution of 



equation fl44[ ), which coincides with the well known equation for the magnetic field of an 
isolated vortex in ordinary superconductors. We can use now these solutions and calculate 
the energy of an isolated quark vortex in the form 




£f (47) 



As we see from this expression the energy of quark vortex depends logarithmically on the 
renormalized charge q. 

In concluding this section let us add the following: If the quark matter occupies a 
spherical volume of radius a, and the vortex passes through the center of this volume then 
the force lines of the magnetic vortex outside of this volume appear to have the dipole form, 
whereas the force lines for the gluomagnetic field have to be confined within this volume 
because of the confinement condition G 8 (a) = 0. Outside of the quark matter region we 
have only the ordinary magnetic field. We also note that we can use these results in the 
calculation of the distribution of magnetic and gluomagnetic fields in a neutron star if its 
core consists of superconducting quark matter. 



IV. CONCLUSION 

In this article we derived the Ginzburg-Landau equations with account for rotated elec- 
tromagnetism. From the form of these equations written for the physical fields B and B$ it 
is obvious that they are coupled, i.e. the appearance of one of them leads to the immediate 
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appearance of the other. We have solved these equations for an isolated quantum vortex in 
quark matter and have shown that the centers of electromagnetic and gluomagnetic vortices 
coincide. The distribution around the continuation of the center of the vortex in quark mat- 

— * — * 

ter for the fields B and B 8 has the same form as for ordinary superconductors. Nevertheless, 
the flux of the magnetic vortex is 6 times the magnetic flux quantum $ and the penetration 
depth in quark matter is 20 times less than that in hadronic matter. Only the force lines 
of the magnetic field can leave the quark matter volume whereas the gluomagnetic field is 
confined therein. As we have demonstrated in this paper for the case of an isolated vortex 
in quark matter the occurence of electric and color Meissner currents is a consequence of 
the color-superconducting state of quark matter. 
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